Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation.

BACKGROUND: Regulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging. RESULTS: Here we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes. We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters. CONCLUSION: Based on these associations of (linearly) distal genes' transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

Trans-eQTL mapping in gene sets identifies network effects of genetic variants.

Nearly all trait-associated variants identified in genome-wide association studies (GWASs) are noncoding. The cis regulatory effects of these variants have been extensively characterized, but how they affect gene regulation in trans has been the subject of fewer studies because of the difficulty in detecting trans-expression quantitative loci (eQTLs). We developed trans-PCO for detecting trans effects of genetic variants on gene networks. Our simulations demonstrate that trans-PCO substantially outperforms existing trans-eQTL mapping methods. We applied trans-PCO to two gene expression datasets from whole blood, DGN (N = 913) and eQTLGen (N = 31,684), and identified 14,985 high-quality trans-eSNP-module pairs associated with 197 co-expression gene modules and biological processes. We performed colocalization analyses between GWAS loci of 46 complex traits and the trans-eQTLs. We demonstrated that the identified trans effects can help us understand how trait-associated variants affect gene regulatory networks and biological pathways.

Understanding operon architecture using LEGO bricks.

Many undergraduates struggle to interpret abstract concepts in molecular biology. Modeling can facilitate learning by making these abstract concepts tangible. Here, we present an exercise based on the lac operon designed for undergraduate students using LEGO bricks. The lac operon is a classic example of transcriptional regulation taught in a variety of undergraduate biology courses and is fundamental to understanding the regulation of gene expression. This easy-to-implement active learning exercise demonstrates how the various components of the lac operon are oriented under a variety of nutritional conditions to control gene expression. In addition, higher-order concepts, such as the effect of mutation on lac operon expression, can be readily modeled. Overall, students not only found this exercise to be enjoyable but also helpful as a tool to engage with this course material.

Cis-regulatory control of transcriptional timing and noise in response to estrogen.

Cis-regulatory elements control transcription levels, temporal dynamics, and cell-cell variation or transcriptional noise. However, the combination of regulatory features that control these different attributes is not fully understood. Here, we used single-cell RNA-seq during an estrogen treatment time course and machine learning to identify predictors of expression timing and noise. We found that genes with multiple active enhancers exhibit faster temporal responses. We verified this finding by showing that manipulation of enhancer activity changes the temporal response of estrogen target genes. Analysis of transcriptional noise uncovered a relationship between promoter and enhancer activity, with active promoters associated with low noise and active enhancers linked to high noise. Finally, we observed that co-expression across single cells is an emergent property associated with chromatin looping, timing, and noise. Overall, our results indicate a fundamental tradeoff between a gene's ability to quickly respond to incoming signals and maintain low variation across cells.

Salicylic Acid Spray Delays Sand Pear Fruit Senescence during Room Temperature Shelf Life by Regulating Antioxidant Capacity and Senescence-Related Genes.

'Whangkeumbae' (Pyrus pyrifolia) is a variety of sand pear fruit well-known for its smooth surface and good taste. However, the fruit quality is adversely affected by postharvest ethylene production. Therefore, improving postharvest shelf life by regulating fruit senescence is critical to promoting the 'Whangkeumbae' fruit industry. Here, we investigated the effect of salicylic acid (SA) spray on fruit senescence in sand pears during room temperature shelf life. Exogenous SA reduced polyphenol oxidase (PPO) activity and malondialdehyde (MDA) content during room temperature shelf life. Additionally, SA effectively maintained the fruit skin coloration and increased the activity of antioxidant enzymes, such as superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), and ascorbate peroxidase (APX). SA treatment inhibited PpPPO1 expression and upregulated PpSOD1, PpAPX6, and PpGST2 expression. Furthermore, SA application downregulated the expression of PpACO2, PpEIN3a, PpNCED1, and PpAOC2, while upregulating PpNPR-1, PpTAR2, and PpCOMT1 during room temperature shelf life. SA treatment also influenced cell wall metabolism and modification genes by inhibiting PpPG1, PpPME2, and PpCEL3 and inducing PpPGIP1 expression. Additionally, SA treatment affected sugar and acid metabolism genes and increased the expression of PpSPS1, PpSUS1, PpSOT1, PpTMT4, PpSWEET15, and PpcyNAD-MDH, but suppressed the expression of PpcyNADP-ME. The Pearson correlation analysis indicated that PPO activity and MDA content were positively correlated with the expression of PpPPO1, PpACO2, PpEIN3a, PpNCED1, PpAOC2, PpPG1, PpPME2, PpCEL3, and PpcyNDA-MDH. Conversely, these factors were negatively associated with the activities of SOD, POD, CAT, and APX, as well as the expression levels of PpSOD1, PpPOD1, PpCAT1, PpAPX6, PpGST2, PpNPR-1, PpTAR2, PpCOMT1, PpPGIP1, PpSPS1, PpSUS1, PpSOT1, PpTMT4, PpSWEET15, and PpcyNAD-MDH. Our results reveal that exogenous SA could delay fruit senescence in sand pear fruit by regulating various biochemical and molecular mechanisms and can be used to effectively extend fruit shelf life during room temperature storage. However, further research is necessary to determine whether the fruits sprayed with SA are suitable for direct human consumption.

Two-way feedback between chromatin compaction and histone modification state explains Saccharomyces cerevisiae heterochromatin bistability.

Compact chromatin is closely linked with gene silencing in part by sterically masking access to promoters, inhibiting transcription factor binding and preventing polymerase from efficiently transcribing a gene. However, a broader hypothesis suggests that chromatin compaction can be both a cause and a consequence of the locus histone modification state, with a tight bidirectional interaction underpinning bistable transcriptional states. To rigorously test this hypothesis, we developed a mathematical model for the dynamics of the HMR locus in Saccharomyces cerevisiae, that incorporates activating histone modifications, silencing proteins, and a dynamic, acetylation-dependent, three-dimensional locus size. Chromatin compaction enhances silencer protein binding, which in turn feeds back to remove activating histone modifications, leading to further compaction. The bistable output of the model was in good agreement with prior quantitative data, including switching rates from expressed to silent states (and vice versa), and protein binding/histone modification levels within the locus. We then tested the model by predicting changes in switching rates as the genetic length of the locus was increased, which were then experimentally verified. Such bidirectional feedback between chromatin compaction and the histone modification state may be a widespread and important regulatory mechanism given the hallmarks of many heterochromatic regions: physical chromatin compaction and dimerizing (or multivalent) silencing proteins.

Shared functions of Fe-S cluster assembly and Moco biosynthesis.

Molybdenum cofactor (Moco) biosynthesis is a complex process that involves the coordinated function of several proteins. In the recent years it has become evident that the availability of Fe-S clusters play an important role for the biosynthesis of Moco. First, the MoaA protein binds two [4Fe-4S] clusters per monomer. Second, the expression of the moaABCDE and moeAB operons is regulated by FNR, which senses the availability of oxygen via a functional [4Fe-4S] cluster. Finally, the conversion of cyclic pyranopterin monophosphate to molybdopterin requires the availability of the L-cysteine desulfurase IscS, which is an enzyme involved in the transfer of sulfur to various acceptor proteins with a main role in the assembly of Fe-S clusters. In this review, we dissect the dependence of the production of active molybdoenzymes in detail, starting from the regulation of gene expression and further explaining sulfur delivery and Fe-S cluster insertion into target enzymes. Further, Fe-S cluster assembly is also linked to iron availability. While the abundance of selected molybdoenzymes is largely decreased under iron-limiting conditions, we explain that the expression of the genes is dependent on an active FNR protein. FNR is a very important transcription factor that represents the master-switch for the expression of target genes in response to anaerobiosis. Moco biosynthesis is further directly dependent on the presence of ArcA and also on an active Fur protein.

ASCL1 promotes Scrt2 expression in the neural tube.

ASCL1 is a transcription factor that directs neural progenitors towards lineage differentiation. Although many of the molecular mechanisms underlying its action have been described, several of its targets remain unidentified. We identified in the chick genome a putative enhancer (cE1) upstream of the transcription factor Scratch2 (Scrt2) locus with a predicted heterodimerization motif for ASCL1 and POU3F2. In this study, we investigated the role of ASCL1 and this enhancer in regulating the expression of the Scrt2 in the embryonic spinal cord. We confirmed that cE1 region interacted with the Scrt2 promoter. cE1 was sufficient to mediate ASCL1-driven expression in the neural tube through the heterodimerization sites. Moreover, Scrt2 expression was inhibited when we removed cE1 from the genome. These findings strongly indicate that ASCL1 regulates Scrt2 transcription in the neural tube through cE1.

Protein degradation of Lsd1 is mediated by Bre1 yet opposed by Lsd1-interacting lncRNAs during fly follicle development.

Tissue development, homeostasis, and repair all require efficient progenitor expansion. Lysine-specific demethylase 1 (Lsd1) maintains plastic epigenetic states to promote progenitor proliferation while overexpressed Lsd1 protein causes oncogenic gene expression in cancer cells. However, the precise regulation of Lsd1 protein expression at the molecular level to drive progenitor differentiation remains unclear. Here, using Drosophila melanogaster oogenesis as our experimental system, we discovered molecular machineries that modify Lsd1 protein stability in vivo. Through genetic and biochemical analyses, an E3 ubiquitin ligase, Bre1, was identified as required for follicle progenitor differentiation, likely by mediating Lsd1 protein degradation. Interestingly, specific Lsd1-interacting long non-coding RNAs (LINRs) were found to antagonize Bre1-mediated Lsd1 protein degradation. The intricate interplay discovered among the Lsd1 complex, LINRs and Bre1 provides insight into how Lsd1 protein stability is fine-tuned to underlie progenitor differentiation in vivo.

Functions of RNAi Pathways in Ribosomal RNA Regulation.

Argonaute proteins, guided by small RNAs, play crucial roles in gene regulation and genome protection through RNA interference (RNAi)-related mechanisms. Ribosomal RNAs (rRNAs), encoded by repeated rDNA units, constitute the core of the ribosome being the most abundant cellular transcripts. rDNA clusters also serve as sources of small RNAs, which are loaded into Argonaute proteins and are able to regulate rDNA itself or affect other gene targets. In this review, we consider the impact of small RNA pathways, specifically siRNAs and piRNAs, on rRNA gene regulation. Data from diverse eukaryotic organisms suggest the potential involvement of small RNAs in various molecular processes related to the rDNA transcription and rRNA fate. Endogenous siRNAs are integral to the chromatin-based silencing of rDNA loci in plants and have been shown to repress rDNA transcription in animals. Small RNAs also play a role in maintaining the integrity of rDNA clusters and may function in the cellular response to rDNA damage. Studies on the impact of RNAi and small RNAs on rRNA provide vast opportunities for future exploration.

Timeless noncoding DNA contains cell-type preferential enhancers important for proper Drosophila circadian regulation.

To address the contribution of transcriptional regulation to Drosophila clock gene expression and to behavior, we generated a series of CRISPR-mediated deletions within two regions of the circadian gene timeless (tim), an intronic E-box region and an upstream E-box region that are both recognized by the key transcription factor Clock (Clk) and its heterodimeric partner Cycle. The upstream deletions but not an intronic deletion dramatically impact tim expression in fly heads; the biggest upstream deletion reduces peak RNA levels and tim RNA cycling amplitude to about 15% of normal, and there are similar effects on tim protein (TIM). The cycling amplitude of other clock genes is also strongly reduced, in these cases due to increases in trough levels. These data underscore the important contribution of the upstream E-box enhancer region to tim expression and of TIM to clock gene transcriptional repression in fly heads. Surprisingly, tim expression in clock neurons is only modestly affected by the biggest upstream deletion and is similarly affected by a deletion of the intronic E-box region. This distinction between clock neurons and glia is paralleled by a dramatically enhanced accessibility of the intronic enhancer region within clock neurons. This distinctive feature of tim chromatin was revealed by ATAC-seq (assay for transposase-accessible chromatin with sequencing) assays of purified neurons and glia as well as of fly heads. The enhanced cell type-specific accessibility of the intronic enhancer region explains the resilience of clock neuron tim expression and circadian behavior to deletion of the otherwise more prominent upstream tim E-box region.

Yeast Crf1p is an activator with different roles in regulation of target genes.

Under stress conditions, ribosome biogenesis is downregulated. This process requires that expression of ribosomal RNA, ribosomal protein, and ribosome biogenesis genes be controlled in a coordinated fashion. The mechanistic Target of Rapamycin Complex 1 (mTORC1) participates in sensing unfavorable conditions to effect the requisite change in gene expression. In Saccharomyces cerevisiae, downregulation of ribosomal protein genes involves dissociation of the activator Ifh1p in a process that depends on Utp22p, a protein that also functions in pre-rRNA processing. Ifh1p has a paralog, Crf1p, which was implicated in communicating mTORC1 inhibition and hence was perceived as a repressor. We focus here on two ribosomal biogenesis genes, encoding Utp22p and the high mobility group protein Hmo1p, both of which are required for communication of mTORC1 inhibition to target genes. Crf1p functions as an activator on these genes as evidenced by reduced mRNA abundance and RNA polymerase II occupancy in a crf1Δ strain. Inhibition of mTORC1 has distinct effects on expression of HMO1 and UTP22; for example, on UTP22, but not on HMO1, the presence of Crf1p promotes the stable depletion of Ifh1p. Our data suggest that Crf1p functions as a weak activator, and that it may be required to prevent re-binding of Ifh1p to some gene promoters after mTORC1 inhibition in situations when Ifh1p is available. We propose that the inclusion of genes encoding proteins required for mTORC1-mediated downregulation of ribosomal protein genes in the same regulatory circuit as the ribosomal protein genes serves to optimize transcriptional responses during mTORC1 inhibition.

Epigenetic dissection of human blood group genes reveals regulatory elements and detailed characteristics of KEL and four other loci.

BACKGROUND: Blood typing is essential for safe transfusions and is performed serologically or genetically. Genotyping predominantly focuses on coding regions, but non-coding variants may affect gene regulation, as demonstrated in the ABO, FY and XG systems. To uncover regulatory loci, we expanded a recently developed bioinformatics pipeline for discovery of non-coding variants by including additional epigenetic datasets. METHODS: Multiple datasets including ChIP-seq with erythroid transcription factors (TFs), histone modifications (H3K27ac, H3K4me1), and chromatin accessibility (ATAC-seq) were analyzed. Candidate regulatory regions were investigated for activity (luciferase assays) and TF binding (electrophoretic mobility shift assay, EMSA, and mass spectrometry, MS). RESULTS: In total, 814 potential regulatory sites in 47 blood-group-related genes were identified where one or more erythroid TFs bound. Enhancer candidates in CR1, EMP3, ABCB6, and ABCC4 indicated by ATAC-seq, histone markers, and co-occupancy of 4 TFs (GATA1/KLF1/RUNX1/NFE2) were investigated but only CR1 and ABCC4 showed increased transcription. Co-occupancy of GATA1 and KLF1 was observed in the KEL promoter, previously reported to contain GATA1 and Sp1 sites. TF binding energy scores decreased when three naturally occurring variants were introduced into GATA1 and KLF1 motifs. Two of three GATA1 sites and the KLF1 site were confirmed functionally. EMSA and MS demonstrated increased GATA1 and KLF1 binding to the wild-type compared to variant motifs. DISCUSSION: This combined bioinformatics and experimental approach revealed multiple candidate regulatory regions and predicted TF co-occupancy sites. The KEL promoter was characterized in detail, indicating that two adjacent GATA1 and KLF1 motifs are most crucial for transcription.

Major transcription factor families at the nexus of regulating abiotic stress response in millets: a comprehensive review.

Millets stand out as a sustainable crop with the potential to address the issues of food insecurity and malnutrition. These small-seeded, drought-resistant cereals have adapted to survive a broad spectrum of abiotic stresses. Researchers are keen on unravelling the regulatory mechanisms that empower millets to withstand environmental adversities. The aim is to leverage these identified genetic determinants from millets for enhancing the stress tolerance of major cereal crops through genetic engineering or breeding. This review sheds light on transcription factors (TFs) that govern diverse abiotic stress responses and play role in conferring tolerance to various abiotic stresses in millets. Specifically, the molecular functions and expression patterns of investigated TFs from various families, including bHLH, bZIP, DREB, HSF, MYB, NAC, NF-Y and WRKY, are comprehensively discussed. It also explores the potential of TFs in developing stress-tolerant crops, presenting a comprehensive discussion on diverse strategies for their integration.

Gene expression regulation by modulating Hfq expression in coordination with tailor-made sRNA-based knockdown in Escherichia coli.

The tailor-made synthetic sRNA-based gene expression knockdown system has demonstrated its efficacy in achieving pathway balancing in microbes, facilitating precise target gene repression and fine-tuned control of gene expression. This system operates under a competitive mode of gene regulation, wherein the tailor-made synthetic sRNA shares the intrinsic intracellular Hfq protein with other RNAs. The limited intracellular Hfq amount has the potential to become a constraining factor in the post-transcription regulation of sRNAs. To enhance the efficiency of the tailor-made sRNA gene expression regulation platform, we introduced an Hfq expression level modulation-coordinated sRNA-based gene knockdown system. This system comprises tailor-made sRNA expression cassettes that produce varying Hfq expression levels using different strength promoters. Modulating the expression levels of Hfq significantly improved the repressing capacity of sRNA, as evidenced by evaluations with four fluorescence proteins. In order to validate the practical application of this system, we applied the Hfq-modulated sRNA-based gene knockdown cassette to Escherichia coli strains producing 5-aminolevulinic acid and L-tyrosine. Diversifying the expression levels of metabolic enzymes through this cassette resulted in substantial increases of 74.6% in 5-aminolevulinic acid and 144% in L-tyrosine production. Tailor-made synthetic sRNA-based gene expression knockdown system, coupled with Hfq copy modulation, exhibits potential for optimizing metabolic fluxes through biosynthetic pathways, thereby enhancing the production yields of bioproducts.

Why an integrated view of gene expression studies on hematopoiesis in mouse aging is better than the sum of their parts.

Globally, the human population is aging, with an increased proportion of people in "old age" (over 60 years). This trend leads to a growing demand in aging research, stimulating studies in animal models such as mice, fish, and invertebrates. Recently, we published a research summary on the aging of hematopoietic stem cells (HSCs) in C57BL/6 mice based on 12 gene expression datasets. Here, I discuss in greater detail the added value of taking an integrated view, rather than considering each publication separately, to determine genes involved in aging. Considerable variation exists between lists of differentially expressed (DE) genes in HSCs, comparing young and old mice. This variation can result from factors such as inconsistent definitions of "young" and "old", technical variations and variations between laboratory mouse strains. We previously demonstrated that the variation between gene lists could be circumvented by forming a unified list of DE genes-the "aging list"-with citation indexes attached. The most frequently detected DE genes [approximately 200 most cited, which we named the "aging signature" (AS)] were highly consistent across publications. Gene Ontology classification of the AS list identified additional sources of variation between studies: one comes from the specifics of how the data are collected and analyzed; another comes from inconsistencies between how we define the gene categories. As discussed, overcoming these variations is the next challenge toward an integral approach to our systematic knowledge of the aging process.

Enricherator: A Bayesian Method for Inferring Regularized Genome-wide Enrichments from Sequencing Count Data.

A pervasive question in biological research studying gene regulation, chromatin structure, or genomics is where, and to what extent, does a signal of interest arise genome-wide? This question is addressed using a variety of methods relying on high-throughput sequencing data as their final output, including ChIP-seq for protein-DNA interactions,1 GapR-seq for measuring supercoiling,2 and HBD-seq or DRIP-seq for R-loop positioning.3,4 Current computational methods to calculate genome-wide enrichment of the signal of interest usually do not properly handle the count-based nature of sequencing data, they often do not make use of the local correlation structure of sequencing data, and they do not apply any regularization of enrichment estimates. This can result in unrealistic estimates of the true underlying biological enrichment of interest, unrealistically low estimates of confidence in point estimates of enrichment (or no estimates of confidence at all), unrealistic gyrations in enrichment estimates at very close (<10 bp) genomic loci due to noise inherent in sequencing data, and in a multiple-hypothesis testing problem during interpretation of genome-wide enrichment estimates. We developed a tool called Enricherator to infer genome-wide enrichments from sequencing count data. Enricherator uses the variational Bayes algorithm to fit a generalized linear model to sequencing count data and to sample from the approximate posterior distribution of enrichment estimates (https://github.com/jwschroeder3/enricherator). Enrichments inferred by Enricherator more precisely identify known binding sites in cases where low coverage between binding sites leads to false-positive peak calls in these noisy regions of the genome; these benefits extend to published datasets.

Crosstalk interactions between transcription factors ERRα and PPARα assist PPARα-mediated gene expression.

OBJECTIVE: The peroxisome proliferator-activated receptor α (PPARα) is a transcription factor driving target genes involved in fatty acid ß-oxidation. To what extent various PPARα interacting proteins may assist its function as a transcription factor is incompletely understood. An ORFeome-wide unbiased mammalian protein-protein interaction trap (MAPPIT) using PPARα as bait revealed a PPARα-ligand-dependent interaction with the orphan nuclear receptor estrogen-related receptor α (ERRα). The goal of this study was to characterize the nature of the interaction in depth and to explore whether it was of physiological relevance. METHODS: We used orthogonal protein-protein interaction assays and pharmacological inhibitors of ERRα in various systems to confirm a functional interaction and study the impact of crosstalk mechanisms. To characterize the interaction surfaces and contact points we applied a random mutagenesis screen and structural overlays. We pinpointed the extent of reciprocal ligand effects of both nuclear receptors via coregulator peptide recruitment assays. On PPARα targets revealed from a genome-wide transcriptome analysis, we performed an ERRα chromatin immunoprecipitation analysis on both fast and fed mouse livers. RESULTS: Random mutagenesis scanning of PPARα's ligand-binding domain and coregulator profiling experiments supported the involvement of (a) bridging coregulator(s), while recapitulation of the interaction in vitro indicated the possibility of a trimeric interaction with RXRα. The PPARα·ERRα interaction depends on 3 C-terminal residues within helix 12 of ERRα and is strengthened by both PGC1α and serum deprivation. Pharmacological inhibition of ERRα decreased the interaction of ERRα to ligand-activated PPARα and revealed a transcriptome in line with enhanced mRNA expression of prototypical PPARα target genes, suggesting a role for ERRα as a transcriptional repressor. Strikingly, on other PPARα targets, including the isolated PDK4 enhancer, ERRα behaved oppositely. Chromatin immunoprecipitation analyses demonstrate a PPARα ligand-dependent ERRα recruitment onto chromatin at PPARα-binding regions, which is lost following ERRα inhibition in fed mouse livers. CONCLUSIONS: Our data support the coexistence of multiple layers of transcriptional crosstalk mechanisms between PPARα and ERRα, which may serve to finetune the activity of PPARα as a nutrient-sensing transcription factor.

Function and regulation of plant ARGONAUTE proteins in response to environmental challenges: a review.

Environmental stresses diversely affect multiple processes related to the growth, development, and yield of many crops worldwide. In response, plants have developed numerous sophisticated defense mechanisms at the cellular and subcellular levels to react and adapt to biotic and abiotic stressors. RNA silencing, which is an innate immune mechanism, mediates sequence-specific gene expression regulation in higher eukaryotes. ARGONAUTE (AGO) proteins are essential components of the RNA-induced silencing complex (RISC). They bind to small noncoding RNAs (sRNAs) and target complementary RNAs, causing translational repression or triggering endonucleolytic cleavage pathways. In this review, we aim to illustrate the recently published molecular functions, regulatory mechanisms, and biological roles of AGO family proteins in model plants and cash crops, especially in the defense against diverse biotic and abiotic stresses, which could be helpful in crop improvement and stress tolerance in various plants.

Promoter regulatory mode evolution enhances the high multidrug resistance of tmexCD1-toprJ1.

Antibiotic resistance could rapidly emerge from acquiring the mobile antibiotic resistance genes, which are commonly evolved from an intrinsic gene. The emergence of the plasmid-borne mobilized efflux pump gene cluster tmexCD1-toprJ1 renders the last-resort antibiotic tigecycline ineffective, although its evolutionary mechanism remains unclear. In this study, we investigate the regulatory mechanisms of the progenitor NfxB-MexCD-OprJ, a chromosomally encoded operon that does not mediate antibiotic resistance in the wild-type version, and its homologs, TNfxB1-TMexCD1-TOprJ1 mediating high-level tigecycline resistance, and TNfxB3-TMexCD3-TOprJ1. Mechanistic studies demonstrated that in nfxB-mexCD-oprJ, MexCD expression was under a weaker promoter, PmexC and inhibited by a strong repressor NfxB. For tmexCD1-toprJ1, TMexCD1 was highly expressed owing to the presence of a strong promoter, PtmexC1, and an inactive suppressor, TNfxB1, with a T39R mutation that rendered it unable to bind to promoter DNA. In tnfxB3-tmexCD3-toprJ1b, TMexCD3 expression was intermediate because of the local regulator TNfxB3, which binds to two inverted repeat sequences of PtmexC. Additionally, TNfxB3 exhibited lower protein expression and weaker DNA binding affinity than its ancestor NfxB, together with their promoter activities difference explaining the different expression levels of tmexCD-toprJ homologs. Distinct fitness burdens on these homologs-carrying bacteria were observed due to the corresponding expression level, which might be associated with their global prevalence. In summary, our data depict the mechanisms underlying the evolution and dissemination of an important mobile antibiotic resistance gene from an intrinsic chromosomal gene.IMPORTANCEAs antibiotic resistance seriously challenges global health, tigecycline is one of the few effective drugs in the pipeline against infections caused by multidrug-resistant pathogens. Our previous work identified a novel tigecycline resistance efflux pump gene cluster tmexCD1-toprJ1 in animals and humans, together with its various variants, a rising clinical concern. Herein, this study focused on how the local regulation modes of tmexCD1-toprJ1 evolved to a highly expressed efflux pump. Through comparative analysis between three tnfxB-tmexCD-toprJ homologs and their progenitor nfxB-mexCD-oprJ, modes, we demonstrated the evolutionary dynamics from a chromosomal silent gene to an active state. We found the de-repression of the local regulator and an increase of the promoter activity work together to promote a high production of drug efflux machines and enhance multidrug resistance. Our findings revealed that TMexCD1-TOprJ1 adopts a distinct evolutionary path to achieve higher multidrug resistance, urgently needing tight surveillance.